ENERGY TRANSFER IN A FAST-SLOW HAMILTONIAN 

SYSTEM 



DMITRY DOLGOPYAT AND CARLANGELO LIVERANI 

Abstract. We consider a finite region of a lattice of weakly interacting geo- 
desic flows on manifolds of negative curvature and we show that, when rescaling 
the interactions and the time appropriately, the energies of the flows evolve 
according to a non linear diffusion equation. This is a first step toward the 
derivation of macroscopic equations from a Hamiltonian microscopic dynamics 
in the case of weakly coupled systems. 

1. Introduction 

One of the central problems in the study of non-equilibrium statistical physics is 
the derivation of transport equations for conserved quantities, in particular energy 
transport, from first principles, (see [6], and references therein, or [30], for a more 
general discussion on the derivation of macroscopic equations from microscopic 
dynamics). 

Lately several results have appeared trying to bring new perspective to the above 
problem in a collective effort to attack the problem from different points of views. 
Let us just mention, as examples, papers considering stochastic models [3, 4, 5], 
approaches starting from kinetic equations or assuming extra hypotheses [2, 26, 7] 
or papers trying to take advantage of the point of view and results developed in 
the field of Dynamical Systems [16, 13, 14, 15, 8, 9, 29]. This paper belongs to 
the latter category but it is closely related to results obtained for stochastic models 
(e.g., [25]). _ 

We consider a microscopic dynamics determined by a (classical) Hamiltonian 
describing a finite number of weakly interacting strongly chaotic systems and we 
explore the following strategy to derive a macroscopic evolution: first one looks at 
times for which we have an effective energy exchange between interacting systems, 
then takes the limit for the strength of the interaction going to zero and hopes to 
obtain a self-contained equation describing the evolution of the energies only. We 
call such an equation mesoscopic since most of the degrees of freedom have been av- 
eraged out. Second, one performs on such a mesoscopic equation a thermodynamic 
limit to obtain a macroscopic evolution. In particular, one can consider a scaling 
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limit of the diffusive type in order to obtain a non linear heat equation as in the 
case of the so called hydrodynamics limit for particle systems, see [22, 31] for more 
details. A similar strategy has been carried out, at a heuristic level, in [19, 20]. 

The first step of such a program is accomplished in this paper. It is interesting 
to note that the mesoscopic equation that we obtain seems to have some very nat- 
ural and universal structure since it holds also when starting from different models. 
Indeed, essentially the same equation is obtained in [25] for a system of coupled 
nonlinear oscillators in the presence of an energy preserving randomness. In addi- 
tion, such an equation is almost identical to the one studied in [31] apart from the 
necessary difference that the diffusion is a degenerate one. Indeed, since it describes 
the evolution of energies, and energies are positive, the diffusion coefficients must 
necessarily be zero when one energy is zero. 

Since, due to the weak interaction, the energies vary very slowly, once the time 
is rescalcd so that the energies evolve on times of order one all the other variables 
will evolve extremely fast. Thus our result is an example of averaging theory for 
slow-fast systems. Yet, in our case the currents have zero average which means 
that standard averaging theory (such as, e.g. [18]) cannot suffice. It is necessary to 
look at longer times when the fluctuations play a fundamental role. The study of 
such longer times can in principle be accomplished thanks to the theory developed 
in [12]. 

Unfortunately, the results in [12] do not apply directly and we are forced to a 
roundabout in order to obtain the wanted result. Not surprisingly, the trouble takes 
place at low energies. We have thus to investigate with particular care the behavior 
of the system at low energies. In particular, we prove that the probability for any 
particle to reach zero energy, in the relevant time scale, tends to zero. 

The structure of the paper is as follows: section 2 contains the precise description 
of the microscopic model and the statement of the results. Section 3 describes the 
logic of the proof at a non technical level and points out the technical difficulties 
that must be overcame to make the argument rigorous. In the following section 
we show how to modify the dynamics at low energies in such a way that existing 
results can be applied. Then, in section 5, we investigate the modified dynamics 
and show that its accumulation points satisfy a mesoscopic equation of the wanted 
type. In section 6 we compute explicitly the properties of the coefficients of the 
limit equation for the modified dynamics and in section 7 we use this knowledge to 
show that the equation has a unique solution, hence the modified process converges 
to this solution. In section 8 we discuss the limit equation for the original dynamics 
under the condition that no particle reaches zero in finite time. The fact that this 
condition holds in our model is proven in section 9. The paper ends with two 
appendices. In the first, for the reader convenience, some known results from the 
averaging theory for systems with hyperbolic fast motion are restated in a way 
suitable for our needs. The second appendix contains some boring, but essential, 
computations. 

2. The model and the result 

For d € N, we consider a lattice 1 d and a finite connected region A C Z d . As- 
sociated to each site in A we have the cotangent bundle T*M of a C°° compact 
Riemannian d-dimensional manifold M of strictly negative curvature and the as- 
sociated geodesic flow g l . We have then the phase space M. = (T*M) A and we 
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designate a point as (q x ,p x ), x G A. It is well known that the geodesic flows is a 
Hamiltonian flow. If we define % : T*M — > TM to be the natural isomorphism de- 
fined by w(v) = (i(w), v}g, G being the Riemannian metric, then the Hamiltonian 
reads 1 Hq = X^eA \Px an( i the- symplectic form is given by to = dq A dp. 2 Thus, 
given x € A, the equations of motion take the form (see [27, Section 1] for more 
details) 

q x = i(Px) , 

(2J) ' fit \ 

Px = F(q x ,p x ) , 

where the F is homogeneous in the p x of degree two. Note that, by the Hamiltonian 
structure, e x := \p 2 x is constant in time for each x € A. It is then natural to use the 
variables (q x ,v x ,e x ), where v x := (p x )~^i(p x ) belongs to the unit tangent bundle 
T l M of M. 3 We have then the equations 

Qx V^I^!E 7 

( 2 -2) v x = \/2e^F(q x , v x ) , 

e x = 0, 

where F is homogeneous of second degree in v x . 

Next we want to introduce a small energy exchange between particles. To de- 
scribe such an exchange we introduce a symmetric, non constant, function (poten- 
tial) V <E C°° (M 2 , R) and, for each e > 0, consider the flow g\ determined by the 
Hamiltonian H E = XLeA \p\ + § Tli\ x -y\ = \ V(lx,qy), that is by the equations 

Qx = i(Px) , 

Px = F(q x ,p x ) - e ^2 d qx V(q x ,q y ). 

\y- x \ = l 

Or, alternatively, 4 




/2e^F(q x ,v x ) + —7= V" {v x L x V - V qj: V(q x , q y )} 
(2-3) ^\y^\=X 

e x = -e\[2el 2_j LxV > 
\x-y\=l 

where (VV,w)g = dV(w) and 

(2.4) L x = v x d qx + F(q x ,v x )d Vx 

denotes the generator associated to the geodesic flow of the x particle on T\M. 



i By p% we mean (i{px),i{px)) G(q x ) = {Px,Px) G where G = i»(G). 

2 To be more precise, given the canonical projection 7r(q,p) = g, first define the one form, 
on T(T*M), p j (£) = p(d7r(£)). Then u> := — dui 1 . Given coordinates q on U C M and 
using the coordinates p for the one form p = '^2 i Pidq i £ T*M, we have a; 1 = ^2 i Pidq i and 
10 = Y2i d*li A dp i , as stated. 

"^Clearly e x is the (kinetic) energy of the geodesic flow at x. 

4 In the interacting case one could chose to include the interaction in the energy and define 
e| := ipj + ^ B |=i %)■ This is the choice made in [25]. Yet, in the present context 

\e x — e%\ < |V|oo£, hence the actual choice is irrelevant in the limit e — > and e x turns out to be 
computationally simpler. 
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We will consider random initial conditions of the following type 

E(/(g(0),«(0))= / f(q,v)p(q,v)dm, V/ e C°((T 1 A/) A , R) 

(2.5) J (TiM) A 

e«(0) = ^ > 0, 

where m is the Riemannian measure on (TiM) A and peC 1 . 

Since the currents L X V have zero average with respect to the microcanonical 
measure, one expects that it will take a time of order e -2 in order to see a change 
of energy of order one. It is then natural to introduce the process e x (e~ 2 t) and to 
study the convergence of such a process in the limit e — > 0. 

Our main result is the following. 

Theorem 1. Provided d > 3, the process {e x (e~ 2 t)} defined by (2.3) with initial 
conditions (2.5) converges to a random process {£ x (t)} with values in R A which 
satisfies the stochastic differential equation 

d£ x = ^2 a(£ x ,£ y )dt+ ^ V2(3(£ x ,£ y )dB xy 

(2.6) \x-y\=l \x-y\=l 

£ X (0) = E X >0, 

where B xy are standard Brownian motions which are independent except that 

13 xy By X . 

The coefficients have the following properties: f3 is symmetric and a is antisym- 
metric; /3 e C°([0,oo) 2 ,R+) and[3(a,b) 2 = abG(a,b) where G £ C°°((0, oo) 2 , R+) n 
C 1 ((0,oo) x [0, oo),R + ) and G(a,0) = A(2a)-i for some A> 0. Moreover, 

(2-7) a = {d £i - d £y )/3 2 + ^(S. 1 Sy 1 )^- 

In addition, (2.6) has a unique solution and the probability for one energy to reach 
zero in finite time is zero. 

Remark 2.1. A direct computation shows that the measures with density hp = 
— — l 

n^eA^ 2 e _/3£a! are invariant for the above process for each j3 € R+. Indeed, 
using (2.7), we can write the generator of the process (2.6) in the simple form 

C = W~ E (9 ex -d £y )h /3 2 (d £x -d ey ) 

\x-y\=l 

from which the reversibility of the generator is evident. 

Remark 2.2. The case d = 2 is harder because the second term in (2.7) (which 
otherwise would give the main contribution at small energies) is zero. We believe 
the result to be still true, 5 but a much more detailed (and messy) analysis of (2.6) 
is needed to establish it. As this would considerably increase the length of section 
9 without adding anything really substantial to the paper, we do not pursue such 
matter. 

Remark 2.3. Note that if we could apply [31] to perform the hydrodynamics limit, 
then we would obtain the heat equation. Unfortunately, (2.6) does not satisfies the 
hypotheses of Varadhan's Theorem on several accounts; the most relevant being that 

^That is the fact that zero is unreachalc. 
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the domain where the diffusion takes place is not the all space and a, (3 vanish 
one the boundary of the domain. This is unavoidable as the energy is naturally 
bounded below. Nevertheless, the results of this paper can be considered as a first 
step along the bumpy road to obtaining the heat equation from a purely mechanical 
deterministic model. 6 

Remark 2.4. As a last remark, let us comment on the choice of1 d . This is done 
just to simplify notations: our arguments are of a local nature, hence the structure 
of % d does not play any role in the proof. In particular, one can prove, with exactly 
the same arguments, the following extension of our result. 

Consider a loopless symmetric directed graph G determined by the collection of 
its vertexes V(G) and the collection of its directed edges E(<G). At each vertex v S 
V(&) we associate a mixing geodesic flow as before, consider then the Hamiltonian 

H c= H l p - + l y(o ei ,q e2 ). 

veV(G) (ei,e 2 )eE(G) 

We then have the exact analogous 8 of Theorem 1 for the variables {£v}vev(G) with 
the only difference that the limiting equation now reads 

d£ v = 2J a(£ v ,£ w )dt + ^ V2/3(£ v ,£ w )dB^ V)W ^ 

(2.8) (v,w)£E(G) (v,w)£E(G) 

S v (0) = E v >0, 

where again for each e € E(G), the B e are independent standard Brownian motions 
apart form the fact that Br v>w \ = —Br WtV \. 

An interesting application of the above Remark is the case where G is a complete 
graph (i.e. E(G) — {(vi,V2) ■ vi,V2 € V(G)}) in which case all particles interact 
with each other. 

The rest of the paper is devoted to proving Theorem 1 . Before going in details we 
explain exactly how the various results we are going to derive are collected together 
to prove the Theorem. 

Proof of Theorem 1. Fix T > and let P E be the probability measure, on the 
space C°([0, T], R+), associated to the process {e x {e~ 2 t)} t ^[o,T\ defined by (2.3), 
P 6 .5 to the one defined by (4.1), F$ the one associated to the process {e z ^} with 
z{t) defined by (7.2) and P the one defined by (2.6). Also, let ft 5 = {t 5 > T} 
where t$ = ini{t E M + : mm x ^\£ x (t) < S}. By construction, for each F E C°, 
Ep e (Fln s ) = Ep e ,, (Fin, ) , % (F l n J = E P (F l n J. 

Proposition 7.4 implies that ¥ e .s => Fs and, since is a continuity set for P^, 
lim^oP^tO?) = F s (n s ) = P(ft 4 ). 



"One could object that geodesic motion in negative curvature is not really mechanical. Yet, it 
is possible to construct a bona fide mechanical system which motion is equivalent to a geodesic 
flow in negative curvature [23]. In any case, by Maupertuis' principle, any Hamiltonian system 
can be viewed as a geodesic flow, possibly on a non compact manifold. 

''Directed means that the edges e £ E(G) are ordered pairs (ei,e2), 6 V(G), which is 
interpreted as an edge going from ei to e^. Symmetric means that if (ei,e2) S E(G), then 
(e2,ei) S E(G). Loopless that, for each a £ V(G), (a, a) ^ E(G). This abstract setting reduces 
to the previous one if we choose V(G) = Z d and E(G) = {(x,y) £ Z d X Z d : \x — y\ = 1}. 

^In particular the condition d > 3 refers to the manifolds M not to the lattice or graph. 
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Next, Lemma 9.1, based on estimate (8.1), tells us that lim^o ^M^-) = 0. Thus 

lim lim FJm) = 0. 

Hence P e => P. The informations on the coefficients follow by collecting (8.3), (2.7) 
(proven in Lemma 8.1), Lemmata 6.1 and 6.3. Finally, the uniqueness follows from 
standard results on SDE and the unreachability of zero (Lemma 9.1). □ 

3. Heuristic 

Let us give a sketch of the argument where we ignore all the technical difficulties 
and perform some daring formal computations. 

If we could apply [12, Theorem 7] to equation (2.3) we would obtain a limit- 
ing process characterized by an equation that, after some algebraic manipulations 
detailed in section 7, reads 9 

(3.1) d£ x = a(£ x ,£ y )dt+ V2[3(£ x ,£ y )dB xy 

\x—y\=l \x—y\=l 

where (3(£ x ,£ y ) = /3(£ y ,£ x ) is symmetric and B xy = —B yx are independent stan- 
dard Brownian motions. The marginal of the Gibbs measure on the energy variables 
reads 

dM/3 = ]^[£x e~ p£x d£ x =: hp A x d£ x , 

X 

for each j3 £ [0, oo). Hence we expect such a measure to be invariant for (3.1). 
Even more, on physical grounds (see Lemma 7.1) one expects the process (3.1) to 
be reversible with respect to this measures. A straightforward computation shows 
that the generator associated to the above SDE reads 

£= E a ^ + \ E fiv@e.-de,) 2 , 

|x— y|=l \x—y\=l 

where a xy — a(£ x ,£ y ). (3 xy = j3(£ x ,£ y ). The adjoint with respect to /iq reads 
C* = E { ~ a *v + ^ * 1 O/^v + - d ty )Plv } d C 

\x-y\=l L J 

+ J E Ply{ d £*- d Sv?-T- E fc.(*o««v) 

\x-y\=l \x—y\=l 
\x-y\ = l 

Computing what it means C — C* implies (2.7). 

Remark 3.1. Note that, as expected, a xy = —a yx . Thus d^2 £ x = 0. 

Going to a bit less vague level of analysis, one must notice that since £ x > 0, 
the diffusion equation (3.1) must be degenerate at zero, also it is not clear how 
regular the coefficients a,/3 are. Hence, a priori, it is not even obvious that such 
an equation has a solution and, if so, if such a solution is unique. To investigate 
such an issue it is necessary to obtain some information on the behavior of the 
coefficients at low energies. 



"See Appendix A for a precise statement of the results in [I 2] relevant to our purposes. 
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To this end one can use the explicit formula given in [12, Theorem 7] for the 
diffusion coefficient. This allows to verify that the coefficients are smooth away 
from zero. An explicit, but lengthy, computation yields, for £ x < £ y , 



2 = AS 




(3.2) 



see Lemma 8.1 for details. Thus, in particular, a xy £ x = ^/3 xy + o((3 xy ). 

We will see in section 9 that such a relation, provided d > 2, suffices to prove 
that the set {(£ x ) ■ Tlx — 0} is unreachable and hence to insure that equation 
(3.1) has a unique solution. 

In the rest of the paper we show how to make rigorous the above line of reasoning. 

4. A MODIFIED DYNAMICS 

Since the geodesic flows on manifolds of strictly negative curvature enjoy expo- 
nential decay of correlations [24, 11] we are in a setting very close to the one in 
[12], i.e. we have a slow- fast system in which the fast variables have strong mixing 
properties. 

Unfortunately, the perturbation to the geodesic flows in (2.3) it is not small 
when e x = 0(e), so at low energies one is bound to lose control on the statistical 
properties of the dynamics. The only easy way out would be to prove that the 
limit system spends very little time in configurations in which one particle has low 
energy. 10 If this were the case, then one could first introduce a modified system 
in which one offsets the bad behavior at small energies and then tries to remove 
the cutoff by showing that, in the limit process, the probability to reach very small 
energies is small. We will pursue precisely such a strategy. 

We now define the modified process. Since our equations are Hamiltonian with 
Hamiltonian H = J2 x eAlPx + § J2\ x - y \=i V(q X : %), the simplest approach is to 
modify the kinetic part of the Hamiltonian making it homogeneous of degree one at 
low velocities and decreasing correspondingly the interaction at low energies. More 
precisely, given any two functions tp, <fr <G C°° (K+\{0}, R), consider the Hamiltonians 
H<P,<t> = Y.x£hV( e x) + lY.\x~y\=i ( l ) { e x)4 l {zy)V{ ( lxiqy)i which yield the equations of 
motion 

qx = p \&x 

)i(p x ) + e ^2 4''{ex)(t){e y )V(q x ,q y )i{j>x) 
\x— y|=i 

Px = f'(e x )F(q x ,p x ) + e ^ <f>'(e x )(j)(ey)V(q x ,qy)F(q x ,p x ) 

\x-y\ = l 

-E ^ 0( e ^)0( e a) d q* V 
\y-x\ = l 

e* = -£ ^^{ey) d qx V{i{p x )), 

\y- x \ = l 



10 To investigate low energy situations directly for the coupled geodesic flows seems extremely 
hard: when the kinetic energy is comparable with the potential energy all kind of uncharted 
behaviors, including coexistence of positive entropy and elliptic islands, could occur! 
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with F has in (2.1). 11 Which, in the variables (q x , v x , e x ), reads 
q x = y/2e^ip' (e x )v x + e ^ y/2e^4>'(e x )<p(e y )V(q x , q v )v x 

\x-y\ = l 



(4.1) 



(p\e x )\/2e^F(q x ,v x ) + el ^ (j)\e x )(j)(ey)y n 2e x ~V(q x , q y )F (q x , v x ) 

\x-y\ = l 

\y-x\ = l \y-x\ = l 



e x = - £ 4>{ e x)4>{e-y)V^e^ d qx V{v x ), 

I?/— 3=1=1 

with F as in (2.2). 

Since > 2 = £{v 2 _ i)^) |s _ a|=1 Mfrl ^y^), the manifold - 1 is an 

invariant manifold for the equations (4.1), thus such equations determine a flow in 
the variables (£ x ,e x ) = (q x ,v x ,e x ) £ T X M x K + . 

Finally, we chose ip = ipg and <fi = <fis such that, for all 5 > 0, 

/ x fs if s > <S , , 1 fl if s > S 

4.2 tf S (s) = i — s ; fo(s) = - T -- = s 

where 0,5 is increasing. 

We denote the solution of the above equations (4.1) with initial conditions (£, e) 

by (¥> 5 (t), e e> s (t)). 

Our goal is to apply [12, Theorem 7] to the flow (£ £,<5 (i), e e ' S (t)), see Appendix A 
for a simplified statement (Theorem A.l) adapted to our needs. Before discussing 
the applicability of this a Theorem, there is one last issue we need to take care of: 
the equation for e is clearly degenerate at low energies, this is related to the fact 
that the energies in (4.1) are strictly positive for all times if they are strictly positive 
at time zero. 12 This may create a problem in the limiting process that is bound 
to have a degenerate diffusion coefficient. To handle this problem it turns out to 
be much more convenient to use the variables z x = lne^. In this new variables we 
finally have the equations we are looking for 

q x = u s (z x )v x + | £ (,s{zx)<t>s{e Zy )V(q x ,q y )v x 
\x-y\=l 

v x =<jJs(z x )F(q x ,v x ) + ^- £ Cs(zx)(/)s{e Zv )V(q x ,qy)F(q x ,v x ) 

\x-y\=l 

(4 3) -4 E e-^Me Z *)Me Zy )V q:c V(qx,qy) 

VZ \v-*\=i 

+ 75 v x e-*cf>s(e Zx )<l>6(e Zv )L x V(£ x ,S v ) 
V2 |w-*|=i 

z x = -eV2 e~^Me Zm )&(e z «)L x V(Z x ,ty), 

\y-x\ = l 



^By dq x V we mean the differential of the function V(-, q y ) for any fixed q y . 
12 Indeed, the equation for the energy can be written, near zero, as e x = —ee x G{e^ x ,^), where 
G is a bounded function, hence the solution has the form e x (t) = e x (0)e~ £ fo G ( e ^a=( s ).5( s )) ds . 
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Where L x is as in equation (2.4) and 

'V2ei if z > In 5 



uj s (z) = V2e% (f' s (e z ) = 

(4.4) 

Cs(z) = V2eH' 5 (e z ) -- 



V25 iiz< In 6 - In f 

if z > In S 
4= if z<hi(5-ln8 

'2o 



Remark 4.1. ./Vo£e that we can chose us > v5 and £5 > decreasing. 13 In 
addition, it is possible to arrange that \u>s\c r (i L ,R) < C r e L , where II = (— oo,2L), 
and |C*|c r (K,E) ^ C r 5~i , for each r E N, L,S E R+. VFe lui/Z assume such properties 
in the following. 

Since the total energy is conserved, we can consider equations (4.3) on the set 
(T!M) A x (-oo,L] A for some L > 0. Hence, by the above remark together with 
the (4.2), the vector field in (4.3) has bounded C r norm, as a function of x, z, e, for 
each r G N. 

Let f s (Z,z,e,6) = f' e (l), F B , s (£,z) = (^ 5 (1), z*> 5 (l)), and 
(4.5) Ai(£,z,e) = -y/2 f J2 e ~^M^ T) )M^ iT) )LMQ S (r),e y ' S (r))dr 

■'0. 



>-a/|=i 
then 

(4.6) F c , s ($ t z) = (f s (£,z,e),z + eA s (t,z,e)) . 

Lemma 4.2. Setting Fg(x,z,s) = F e ^(x,z) we have, for each S G (0,1), L > 0, 
F 5 G C°°((T 1 M) A x (-oo,L] A x [0,1]), and ||A*(. I . | e)|| C r ((r i M) A X (_ 00i£] A ) < C r ,s, 
for each r G N, £ G [0,1]. In addition, for each (3 G R+, i/ie probability measure 

d/j,s,e,f> = Z^ x e~ ?Hs ^ + ^ i Zx dqdvdz 

II A> (g,^,z) = ^<^(e z *) + | ^ 0a (e^)^ 2 ")^, g„), 
ieA |o;— i/|=i 

is invariant for F e g. Moreover, for each z G R d and sub-manifold := {z^ = 
the Dynamical System (T,z,Fqj) has a unique SRB measure /ig- 

Proof. The first part of the statement follows from Remark 4.1 and subsequent 
comments together with standard results of existence of solutions and smooth de- 
pendence on the initial data from O.D.E.. The bound on A d is then immediate 
from formula (4.5). 

By the Hamiltonian nature of the equations (4.1) the measures 

dfi s ,p = Z^e-^'^dqdp, 

are invariant for the associated dynamics for each (3 > 0. By changing variables we 
obtain the statement of the Lemma. 

Finally, calling ft the Riemannian measure on T X M we have that [iz = <8>' A '/i is 
a SRB measure for the map £ h-j- / d (£, z,0), which turn out to be the product of 



13 Indeed, 



J min{s,<5} 



t> s (s) = 1 - / p=-(ic. 

/2x 



Remark that once ( $ is chosen all the functions are fixed. 
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the time uis(z x ) maps of the geodesic flow on T X M. The uniqueness of the SRB 
follows by the mixing of the geodesic flows [1] and the fact that the product of 
mixing systems is mixing. □ 

5. Existence of the limit: 6 > 

Wc arc finally ready to consider the limit e — > 0, for the modified dynamics. 

Proposition 5.1. For each S € (0, 1) there exists e$ > such that the Dynamical 
System defined by (4.6) satisfies the hypotheses of Theorem A.l for e G [0, Sg]. 

Hence, the family z £ ' S (e~ 2 t) is tight and its weak accumulation points are a 
solution of the Martingale problem associated to the stochastic differential equation 



(5.1) 

where 



d4 = al(z s )dt + J2°t v (z S )dB y , 
v 

=E / A s x ((f s n^z,0)A%z,Q)d^ z 



(5.2) = 2 / dt 



ne 



/ j z x + z y 

e 2 



\x — w\—l 

\y-w'\=i 

xE(L x V(& s (t),ej(t)) ■ L y V(H y ,Z w ,)) . 
Here E is the expectation with respect to fj, z and Ha^Hco + || (cr 5 ) 2 || c i < oo. 

Proof. First of all notice that the hypotheses on the smoothness of F s j and the 
boundedness of A s are insured by Lemma 4.2. Next, notice that Fo t s(£, z) = 
(fz(£),z) with fz(£) x = g u ' s ^ Z! °'(^x), where g 1 is the geodesic flow on the unit 
tangent bundle T l M, thus the /* are FAE. 14 

Also we have that fi z (A z,0)) = 0. This follows by considering the transfor- 
mation Q(q 7 v) = (q, —v). Indeed &*fi z = fJ-z while, the flow "J^ E associated to (4.3) 
satisfies $j £ o6 = 9o ^jl- O n the other hand, using the antisymmetry of L X V 
with respect to v x , 



A 5 x (Q(0,z,0)=-V2 f V e-*^(e«-)^( C "»)L a! ^o*J i0 oe(0d7 

J i i , 
\x-y\=l 

=V2 f V e-»^(e*-)^(e*»)^Vo^J(0dT 
Jo i_ .., , 



\x-y\=l 

= -4(*- i (e),z,o). 

Thus^(A A '(-,z,0)) = ^ Z (A S (&(■), z,0)) = -n z (A s (^(-),z,0)) = -fi z (A 5 (-, z, 0)), 
by the invariance of the measure. 

The last thing to check is the uniform decay of correlation. Since us > y/d, the 
results in [11, 24] imply 15 that the f z are FAE with uniform exponential decay of 



14 FAEs arc defined in Appendix A. In our case, the abelian action is the one determined by 
the geodesic flows themselves, Xi£\g ti . 

15 [11] proves the exponential decay of correlations for geodesic flows on negatively curved 
surfaces, [ ] extends the results to any negatively curved manifold. 
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correlation. In fact, in Theorem A.l the decay of correlations is meant in a very 
precise technical sense. To see that the results in [24] imply the wanted decay we 
must translate them into the language of standard pairs in which it is formulated 
Theorem A.l. Let us start by stating the result in [24]: let g a be the time a map 
of the geodesic flow on the unit tangent bundle. For each smooth function A let 
||^4||s = |j4||oo + ||9' s 'A||. 0O where d s is the derivative in the weak stable direction. 
Then there exists C, c > such that, for each z and p, A G C 1 , holds true 

(5.3) \E(p-Aog™)-E(A)E(p)\<C\\p\\c4A\\ s e- can . 

Since, setting /*(£) = f s (£, z,Q), /* = Xj,^ 1 ' 2 '', and u}$ is uniformly bounded 
from below, for E(A) = 0, it follows (suppressing, to ease notation, the superscript 
<5) 16 

(5.4) \E(p.Aoft)\<C\A\\\p\\c4A\\ a e- am . 

To see that this is stronger than needed, consider a standard pair I = (D,p). 17 
One can smoothly foliate a e neighborhood of D and define a probability density p e 
supported in it such that 1 1 p e 1 1 c 1 < Ce~ 2 , while || c i < C when p e is restricted to 
a leaf of the foliation. Thanks to the a-H61der regularity and the absolute continuity 
of the weak stable foliation, one can take p £ so that 

\E t (A)-E(p E A)\ <Ce a \\A\\ s . 

Accordingly, 

\E t (A o f:)\ < |E(p E ■Aof:)\ + Ce a \\Aof:\\ s <C {e- 2 e~ can + s a } \\A\\ S 
<Ce~^\\A\y, 

where, in the last equality, we have chosen e = e~ 2 + a . Thus, all the hypotheses of 
Theorem A.l are satisfied and (5.2) follows by a direct computation. □ 

By Theorem A. 1(b), in order to prove that z £ ' S (e~ 2 t) has a limit it suffices to 
prove that (5.1) has a unique solution. This would follow by standard results if 
we knew that a s is locally Lipschitz. In fact, [12] provides also an explicit formula 
for a s . Unfortunately this formula is much more complex than the formula for the 
variance and is quite difficult to investigate. We will avoid a direct computation of 
as and we will instead use the knowledge of the invariant measure to determine it. 
Before doing that a deeper understanding of the variance is required. 

6. Computing the variance 

Let g t be the geodesic flow on the unit cotangent bundle of M. As already noted, 
for each function h, h(£®' (t)) = ho g u ' l '^ Zx ' t {^x) for all a; € A. For convenience let 
us set vj x := uj$(z x ). Also, it turns out to be useful to define two functions of two 
variables: consider two geodesic flows on T 1 M, let (£,77) be the variables of the 
two flows respectively, E the expectation with respect to the Riemannian volume 



Just note that one can write E(p- Ao/™) = E(E(p- Ao f™ \ ^ y ^. x )) and that the relevant norms 
of Pfv^C&O = Piix'iy^x) and A^^^) = A(| x ,/"(^ x )) are bounded by the full norms of p 
and A. One can then apply (5.3) to E(p ■ A o | = E(p f ^ x A 5 ^ x o /^ (Zx) ). Proceeding 

in such a way one variable at a time yields the result. 

1 ^Recall that D is a manifold of fixed size close to the strong unstable one and p a smooth 
density on it. 
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on (T^-M) 2 and L\,L-x the generators associated to the geodesic flow of £ and 77 
respectively, then we define p, p : M 2 — > R by 

/oo 
dt E (L 1 V(g at (^),g bt (r])) ■ la V(£, 77)) , 

p(a, 6) := / dt E (Li 9 M (v)) ■ L 2 V(£, 77)) . 

Also, it is convenient to define 

(6.2) p xy := p(ujs{z x ),ujs{zy)), p xy := p{uss(z x ),Los(z y )). 

Indeed, the understanding of the variance will be reduced shortly to understand- 
ing the properties of p xy . Here is a list of relevant properties whose proof can be 
found in Appendix B. 

Lemma 6.1. The function p is non-positive and C°° for a, b > 0. In addition, 
for each a, b, X > we have p(a,b) = p(b,a) and p(Xa,Xb) = X~ 1 p(a,b). Finally, 
P(a,b) = -fp(a,b). 

Remark 6.2. Note that the previous Lemma implies a 2 p(a,b) = b 2 p(b,a). 
Lemma 6.3. There exists A, B > such that, for all a, b > 0, 

Bab 3 



< 



V> + b 5 ' 



Finally, for all a,b > 0, 

B ab 2 

\9 a p(a, b)\ < ; ad a p(a, b) + bd b p(a, b) = -p(a, b). 

a a + b° 

We are now in the position to derive an helpful formula for the variance. 
Lemma 6.4. The following formula holds true 

-2e-*yct> s (e z *)<t>s(e Zy ) 3 Pxy if\x-y\ = l 

if\x-y\>l. 

Proof. Remembering (5.2), given any two couples of neighboring sites x, w, y,w' 
we want to compute 

dt E (£xn<T**(£x), S""'^)) ■ LyV&vM) ■ 

> 

In fact, remembering the properties of the transformation in the proof of Lemma 
5.1, it suffices to compute the integral on [0, 00). 

Since E^ | q^ x ,v^ x ) = 0, it follows that the above integral is different from 
zero only if x = y or x = w' and w = y. On the other hand if x = y, since g at x g bt 
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is a mixing flow for each a, b > 0, we can write 



m (V(g X) q w )) E (L x y(C XJ &,/)) - w^E <?„,) • L x V(6s> 



/"OG 

= -W— / dtE(L w V( g "" t (£ ai ) t g**' t {Z w ))-L a! V(i x ,Z w )) 



'x Jo 



Jo 

Thus, remembering (4.2), (4.4) and that w x = uis(z x ), 

\x — w\ — l 

and a xy = if |x — y\ > 1. If \x — y\ = 1, then (remembering the symmetry of the 
potential and using Lemma 6.1) 



<j 2 xy = 2<f> 5 (e z *) 2 ct> s (e z vfe- iE ¥ 3L p xy = -2e— »&(e«")&(e*») 3 pxv. 



□ 



7. The limit equation (<$ > 0): structure 

Having gained a good knowledge on the variance we are ready to write the limit 
equation in a more explicit and convenient form. 

We introduce standard Brownian motions B xy indexed by oriented edges, so 
that the motions associated to different non oriented edges are independent and 
B xy = —B yx . Considering the Gaussian processes W x := J2\ x - y \=i fixy(z)B xy we 
have 

(E\ x - w \=iPxUz) 2 1 fovx = y 
E(W x (t)W y (t) \z)=l -p xy {z)fi yx {z) t for \x-y\ = l 

[o for \x - y\ > 1. 

We set 18 

(7.1) (3 xy (z) = V2e-*&(e*")fo(e"») v ^:, 

hence, remembering Lemmata 6.4, 6.1 and equations (6. 2), (4. 4), (4. 2), 

fE|*-«,|=i#L, ifx = y 

{<r S )ly{z) = \ -PxyPyz if IX- 2/1 = 1 

[o ii\x-y\>l. 

Then, we can write (5.1) as 

(7.2) dz s x =a s x (z s )dt + Px V {z S )dB xy . 

\x-y\=l 

Let C be the operator in the Martingale problem associated to the diffusion 
defined by (5.1). 



^This is well defined since p x ^ y > by Lemma 6.1. 
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Lemma 7.1. If the manifold M is d dimensional, then for each /3 > 0, 

zs an invariant measure for the process defined by (7.2). In addition, the process 
(7.2) is reversible. That is, calling the expectation with respect to the above 
invariant measure, 

Ep((p£h) =Ep(h£(p) 
for each smooth real functions tp, h. 

Proof. Recall that Lemma 4.2 gives the invariant measures of the original Dynam- 
ical System. In particular , for each ij> £ C°(R' A ',R) 

\tX5,eA^{z e ' 5 {e-H))) - Lt SM (i/,(z s - s (e-H)))\ < Ce^U 

Thus 

\ii &fi Aip{z e ' 5 {s- 2 t))) - ix Sfi A^ e ' 5 (0)))\ < 2Ce|V|oo. 

Taking the limit e — > along any subsequence leading to an accumulation point 
we see that fJ<s,o,p is an invariant measure for the process (5.1). The claim of the 
Lemma now follows by taking the marginal of Hs,o,i3 in the variables z. 

In the same manner, using the same notation as in the proof of Lemma 5.1, for 
each continuos functions ip,g and converging sequence z ek s(s~k t) we have 

e~ 2 t -e~ 2 t 

Ep(ip{z{t))g(z)) = lira fJ, S ,e k A9 ' i> ° %l,8 ) = , lim Vs,e k A^ ' 9 ) 

= lira Hs,e k A^ ° ' 9 ° © ° KTs) = E p(9 ° e(*(t)ty o 6(z)). 

Since ^ are functions of the z only, it follows ^o0 = ^ : ^o9 = ^ and 

%(^(^(%(z))=E^(^(i))^W). 
Differentiating with respect to t at t = yields the Lemma. □ 

Lemma 7.2. TTie dri/t a s x has the form 

< = E [e-"-^(e"") a ^(e"») 2 Px»] - d Zy [e^Me^M^fp^]} 

|x-J/|=l 

+ ^ E [e" Za; 06(e Za =) 2 ^(e z '') 2 -e- 2 »^(e^)^( e z '') 3 ]^. 

|x-i/|=l 

Proof. The idea to compute the a* is very simple: first compute £ and C* and 
then check what the reversibility condition implies. The operator associated to the 
diffusion (5.1) is given by 

x x,y 
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The adjoint C* with respect to the invariant measures in Lemma 7.1 can then be 
computed by integrating by parts. Setting T x (z) := | — (3(f>s(e Z:c )~ 1 we have 

£ > = - E^ a * + - E a ^ 

X X 

+ \Y, [ d M° S )ly + ^ x d Zy (<T% + T x T y (a s f xy + 8 xy d z T x {a s Y xy ] $ 

xy 

+ E K(° 5 )ly + T y{° S )lv\ « + \ E^sA^M- 

xy xy 

This implies 

y 

and the Lemma follows by direct algebraic computations using Lemma 6.4. □ 

The next result is an obvious fact that is nevertheless of great importance. 

Lemma 7.3. The function % := ^2 X tps(e Zx ) is constant in time. 

Proof. It is useful to notice that, setting ip x := ^ e ^ x ^ and K xy = ip x /3 xy , n xy = n yx . 
By Ito's formula we have 

a"H = E ^xO-xdt + E K xydB xy + ^ E dz ^ x E Ply dt - 

x \x— y|=l x \x-y\ = l 

The second term is zero by the antisymmetry of B xy , thus (using Lemma 7.2 and 
the symmetry of n xy again) 

dn =- E bPxdzxPly ~ ^ydzjxyfiyx] dt 
\x~y\ = l 

+ ^ E [^x 1 -% 1 ]<ydt+\ E PlyW*dt = Q. 
|x-i/|=l |x-j/|=l 

□ 

We conclude with the main result of this section. 

Proposition 7.4. For each 8 > the family z £ - s (e~ 2 t) converges weakly, for e — > 0, 
to i/ie process z(t) determined by the SDE (7.2). 

Proof. From Lemma 7.2 and Lemma 6.1 it follows that a 5 £ C°°, this, together 
with the boundedness and convergence results established in Lemma 5.1 and the 
standard results on the uniqueness of the solution of the SDE, imply that all the 
accumulation points of z € - s (e~ 2 t) must coincide, hence the Proposition. □ 

8. The limit equation (5 = 0): properties and stopping times 

It is natural to consider the stopping time ts := mf{t € R+ : min^gA z x < ln<5}. 
In addition, Lemma 7.3 suggests the convenience of going back to the more physical 
process 8 x {t) = <pj(e*'<* Ar '>) = e z ^ tArs \ 
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Lemma 8.1. For each t < ts, the process £ x satisfies the SDE 

d£ x = Y a(£ x ,£ y )dt + V2j3(£ x ,£ y )dB xy 
y 

where a,/3 € C°°((0, oo) 2 , R) are respectively anti- symmetric and symmetric func- 
tions that satisfy (2.7), (3.2). In addition, ifd>S, then for each constant 

M>nm 1, d_ 2 f' 

(8.1) a(£ x ,£ y )£ x >f3(£ x ,£y) 2 . 

Proof. By Ito's formula and (7.2) we have 19 



(8.2) 



d£~ = 



? z *a x + -e z * % 



3 2 

\x-y\=X 



dt+ ^2 e z * (3 xy dB xy . 

\x-y\=l 



Using (7.1), (6. 2), (4.4) and Lemma 6.1 we can write 
(8.3) 



e Zx A,y = \]2£ xP {^/2£~ x: y/vQ) =: V2f3(£ x ,£ v ). 
Lemma 7.2, equations (4.2), (6.2) and (4.4) yield 



a x = [ d £*P- d £ v P] 

I a: -1/1=1 



d-2 



\x-y\=l 

Using equation (8.2) we finally obtain (2.7) and from Lemma 6.3 follows (3.2). 
Moreover, by Lemma 6.3, 

L -.d a p{V2£Z,^2£~y) 



d£xPxy 
ds y Pxy 



/2£, 
1 

'2£. 

(3(£ x ,£ y ) 2 £. 



l — { P (V2£Z, y/2£y) + ^7£ x d aP {^2£ Xl V^)} 



1£ £ £ 



-z-de x p X y 



Hence 



£ x a(£ x ,£ y ) = (3 2 + £ 2 d £scPxy - £ 2 d £yPxy + ^(3 2 - ( L^£ x £- 1 f3 2 



f3{£ x ,£yf+£ 2 \l + ^\d £i p xv . 



d d — 1 £ x 

The regularity of the coefficients follows from the previous results and some alge- 
braic computations. At last, for £ y > M£ x , 

Bil + M- 1 ) £ x 



£ x a(£ x ,£ y ) > < - 



d d-1 
2 ~ 2AI 



2M (2£ y )h 



19 Here we suppress the <5-dependence since we stop the motion before seeing the region in 
which the dynamics has been modified. 
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On the other hand 

(8.4) (3(£ x ,£ y ) 2 = £ xP {^2£~ x , V^~y) > 
from which the Lemma follows. 



-4 



B\J2£ X 
2£ v 



> 



A£ x 



□ 



9. The limit equation (5 = 0): unreachability of zero energy 

Our last task it to prove that the stopping time t$ tends to infinity when 6 tends 
to zero or, in other words, energy zero is unreachable for the limit equation. 
Fix any T > 0. 

For each subset T C A let us define the energy of the cluster £r := J2 x er^ x - 
Also, for each 6 > 0, n S {1, . . . , |A|}, let us define the stopping times 

t'£ := inf{< e [0,oo) : 3 T C A, |r| = n,£ r (t) < S} AT. 

Note that t$ = t$ A T, where rg is defined at the beginning of section 8. 

Lemma 9.1. Let P be the measure associated to the process (2.6), then 

limPfW <T\) =0. 

Proof. We we will prove that for each 77 > and n <G {1,...|A|} there exists 
$n = finiv)) such that 

P({t£ < T}) < 2- n rj. 

The proof is by (backward) induction. The case n = |A| follows by the energy 
conservation by choosing 8\\\ < 4p. 

Next, suppose the statement true for n + 1 < |A|. It is convenient to define, for 
each r C A the stopped process £ T {t) = ^(tAT^) and the set Q = {r^ > T}. 
Then, for each < S < <5„+i, we have 

P (W <T})<¥ ({t 5 ™ < T} n fi) + 



( 



< 



u U 



rcA 

\\T\=n 



inf £ T (t) < 5 

[0,T] 



\ 



+ 2 -(™+ 1 )^. 



It thus suffices to show that there exists 5„ < 6 n +i such that, for each T C A, 
|r| = n, we have 



inf £ r {t) <S n })<2 



-(|A|+n+l), 



< 



\M\ 2 -(^ 



k te[o,T] 

Let us fix T c A, |r| = n. 

Observe that if fi holds but £ r (t) < |f±i then £ v > > M£ T > M£ x for 

all y ^ r and x E V. In the following we will chose M as in the statement of Lemma 
8.1. 

Next, we define the process Y = \n£p which satisfies 
a(£ x ,£ y )£r - f3(£ x ,£ y ) 2 



2£ 2 



(9.1) dY= Yl 

(x, v )eB(r) 

where B(T) = {(x,y) e A 2 : x € T, y £ T, [ac — 2/| = !}• 



-dt + ^/S^,^)^ 1 ^! 
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Observe that by Corollary 8.1 the drift is positive, indeed 

E (a(£ x ,£ y )£ r - f3(£ x ,£ y ) 2 ) x - (a{£ x ,£ y )£ x - (3(£ Xl £ y ) 2 ) 
2£ 2 -2^ 2£ 2 

(x,y)eB(r) r (x,y)€B(T) r 

In addition, arguing as in (8.4), if £r(t) < |fr^ we have, for some constant C > 0, 



(9.2) S^ 2 0{S x ,Sy) 3 <2- a 



> 

Therefore 



C 

< — 

172 



Y(t A r££ ) > F(0) + / " +1 ^ s/2p(£ x ,£ y )£^dB xy =: M(t). 

' /o (i,Sl)£B(r) 

Note that M is a Martingale. Let r* = inf{* : M(t) < ln<J n+1 } A T. Consider the 
new martingale M.(t) = M(t) — M(t A r*) and the stopping time 

f = inf{f : M(t) < In 5 n - In S n+1 or M(i) > -- lnS n+1 } A T. 

Setting p = P({M(f) = ln<5„}) we obtain 

< p(ln<5„ - ln<S„ + i) - (1 - p)^ln5 n+1 , 



which implies 

(f) = lnU)< 



\n5 n+1 



21n<5„ - hi<5„ + i 

Set 5 n = 6" +1 , a > 1 to be chosen later. The probability that M, starting from 
ln<5„_|_i reaches ln<5„ before reaching | lm5„+i is smaller than (2a — Accord- 
ingly, the probability that the martingale reaches In S n before downcrossing L times 
the interval [ln£„+i, h ln<5 n+ i] is smaller than 1 — (1 — (2a — 1) _1 ) L < a~ x L. On 
the other hand by Doob's inequality the expectation of the number of downcross- 
ing is bounded by 2 _! E((M - ± In <5„ + i) + ). Since M - ±hi(5 n+1 > implies 



£r > v^+T. b y (9-2) follows 

E((M-iln<5„ +1 )+)<C<S;| 1 , 

for some constant C independent on e. From this it immediately follows that the 
probability to have more than L downcrossing is less that L~ 1 6~ +1 . In conclusion, 



Jrf^rW < S n jj < C(a- 1 L + L- 1 5- 1 +l ) 

which yields the wanted estimate by first choosing L 2 = a5~^ x and then setting 
a = C 2 5~l 1 2 2 W+ 2n + A T 1 - 220 □ 

Corollary 9.2. The set {3x : £ x = 0} is inaccessible for the limiting equation. 



2 ^Note that <5 n ~ <5 n+ " +1 for some constant C. So, for large A, <5i is absurdly small. Yet, this 
suffices for our purposes. 
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Appendix A. An averaging Theorem 

In this appendix, for the reader convenience, we recall [12, Theorem 7] stating 
it in reduced generality but in a form directly applicable to our setting. 

Let M be a C°° Riemannian manifold, z £ R d and f z £ Diff°°(M, M) a family 
of partially hyperbolic diffeomorphisms. 21 

We say that {f z } is a family of Anosov elements (FAE) if there exists Abelian 
actions g z ,u t £ M. dc where d c — dimE c , such that f z o g zt = g zt o f z and 
span{5 ti g z , t } = E c . 

Next, we need to discuss decay of correlations that in [12] is meant in a very 
precise technical sense. The basic concept is the one of standard pairs. For the 
present purposes a standard pair can be taken to be a probability measure deter- 
mined by the couple I = (D, p) where D is a C 2 dim(£' u )-dimcnsional manifold D 
close to the strong unstable manifold and a smooth function p £ C 1 (D,M. + ) such 
that J D p = 1- We set Kg(A) = j D Ap. The point is that it is possible to choose 
a set E of manifolds D of uniform bounded diameter and curvature such that, for 
each D £ S, f z D can be covered by a fixed number of elements of E. For each 
C > we consider the set E x = {(£>, p) : DeS, ||p|| c i (D, E) < C} and let E 2 be 
the convex hull of E\ in the space of probability measures. 

It is easy to check that one can chose E and C such that for all I £ E\ there 
exists a family {£i} C E\ such that E| (Aof z ) = c i^ii{A). In addition one can 

insure that any measure with C 1 density with respect to the Riemannian volume 
belongs to the weak closure of E 2 (see [12] for more details). 

We say that the family {f z } has uniform exponential decay of correlations if 
there exists Ci, C 2 > such that, for each z £ M. d there exists probability measures 
p z such that for each n £ N, standard pair £ £ E\ and functions A £ C 1 (M, K.) it 
holds 

\K e (A o f n z ) - p z (A)\ < C ie - c * n \A\ cl . 
Consider now the function F £ C°°{M xM d xl + ,Jl/x M d ), 

(A.l) F(x, z, e) = (f(x, z, s), z + eA(x, z, e)), 

and the associated dynamical systems F e (x,z) = F(x,z,e), such that f(x,z,0) = 
f z (x). Let (x e n (x,z),z e n (x,z)) := F?(x,z). Then for each g £ C r (M,R + ), p(g) = 1 
we can define the measure p g (h) := p(g ■ h) and consider the Dynamical Systems 
(F e ,M x M. d ) with initial conditions z — zq and x distributed according to the 
measure p, g . We can then view z^ as a random variable, clearly E(i/j(z^)) = fJ- g (ip o 
F™), where "ip(x,z) = ifj{z). 

Theorem A.l ([12]). Let F,F e ,f z be defined as in (A.l) and subsequent lines. 
Let f z be FAE with uniform exponential decay of correlation. Suppose that there 
exists £o,C r £ M+ such that sup e<eo ||^4(-, -, < CV and (j, z (A(-, z, 0)) = for 

all z. Also assume that Zq = z* and Xq has a smooth distribution on M as described 
above, then 



By this we mean that, for each fixed z, at each point x £ M the tangent space of T X M 
can be written as E u {x) © E c (x) © E s (x), where the splitting is invariant with respect to the 
dynamics, i.e. d x fE*(x) = E*(f(x)) for * £ {u,c,s}. In addition, there exists constants Ai < 
A2 < A3 < A 4 < A 5 < A 6 , with \2,^ 1 < 1, such that Ai < a(df\ E „) < P/|s s || < A 2 , 
A3 < a(df\ Ec ) < \\df\ Ec \\ < A 4 and A 5 < a(df\ E J < \\df\ E J\ < A 6 , where a(A) = ||>4. 1 1| x . 
22 The integral is with respect the volume form on D induce by the Riemannian metric. 
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a) The family {^f te -2-i} is tight. 

b) There exists functions a 2 E C 1 ^, SL(d, R d )), a 2 > 0, a € C°(R d , R d ) such that 
the accumulation points of {zf te -2-\} are a solution of the Martingale problem 
associated to the diffusion 

dz = adt + adB 
z(0) = z., 

where {Bi}f =1 are independent standard Brownian motions and 

oo „ 

a 2 (z)= V] / A(x,z,0) ® A(fgX,z,0)fj, z (dx). 
'Jm 



n— — oo 

2 



Moreover ||a||co + \\o~ Wc 1 < oo. 

Appendix B. The properties of p xy . 
Here we collect, a bit boring, proofs of the Lemmata concerning p xy . 

Proof of Lemma 6.1. The non-negativity follows from the fact that the quantity 
is an autocorrelation, see footnote 24 for details. By definition 



d^8^p{a, b)= dt t n+m E o g at ® g bt ■ L 2 V) 

J — CO 

/oo 
dt t n+m E ((L r l +1 L™V) ■ L 2 V o g at $ g bt ) . 
-OO 

Applying (5.3) to the above formula yields 

poo 

b)\ < C n , n / dt t n+m e -c m in{aMt < ^ min{a> b} -n-m- 1_ 

Jo 

This proves the smoothness of p. To continue, consider 

/•OO 

/5(Ao, Xb) = / dt E ((Li V) o g aXt g> g bxt ■ L 2 V) = A^pK b) 



by the change of variables t — > Xt. The symmetry follows by a change of variables 
as well. Finally, 



p(a,b)= / dt E ((LiV) o g at ® g bt -L 2 V) 

J —OO 

/OO 
dtE{L x V ■ {L 2 V) o g at ® g bt ) 
-oo 

= b- 1 / dt — E [L X V -Vog at ® g bt ) 
J — oo 

/oo 
dt E (LiV ■ (LiV) o g at ® g 
-oo 



The lemma follows then by the mixing of g at (E> g bt (being the product of two mixing 
flows) and the definition of p. □ 

To continue it is useful to define and study the function T(t) := p(r, 1). 
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Lemma B.l. There exists A,B>0 and fl£l such that 



A 



Bt 



1 + T° ITT 

r'(r) - Dr\ < Bt 2 , 
r'(r) +3At- 4 \ < Bt- 5 , 



Vr > 0, 

Vr G (0, 1] 
Vr > 1. 



Proof. Let us start by assuming r < 1. By setting V(qi) = E(V \ qi,v\), and 
taking care of adding and subtracting that is needed to write convergent integrals, 

/>oo 

r(r) =2 1 dtE (L\V ■ L\V o ff Tt ® gr*)) 

/■OO 

/ fit [E (LiF • LiF o id ® -E((LiF) 2 )] 
'o 

/>oo />rf 

2/ dt I ds [E(L 1 V-LlVog*®g t )-E(L 1 V-LlVog a )] 
Jo Jo 

dt E(LiV ■ LiV o g Tt ) 



=2 



The third term here vanishes since it is the variance of a coboundary. That is, 



dt-EfaV-Vog') 
dt 



dtE(L x V ■ LiV o g Tt ) = t~ 

Also, setting V = V - V, 

/ dt ds [E (L{V ■ L\V og s ® g f ) - E(LiV ■ LfV o g s )\ 
Jo Jo 



0. 



dt E ( L\V ■ LfV o g s 



= o 



ds 



dt = 0{t) 



were we have used (5.3) after conditioning with respect to q\,v\. Thus 

/>oo 

T(r)=2 / dt [E{d qi V -d qi V oid®g l ) ~E{{d qi V) 2 )] 



(B.l) 

+ 0(t) = A + 0(t). 

The fact that A > follows from general theory of mixing flows combined with 
cocycle rigidity of geodesic flows [21, 10]. 24 



2:S 



Here we use the fact that E(i)i (g Vi \ q±, rf) = 1. 
!4 Indeed, for each T > and / 6 C°°, E(/) = 0, 

< E 



(f^fogi-dt \ = 2 ^ dt(T - t)E(/ og t -f) = 2T £ dt E(f o g* ■ /) + O(l). 

Thus the autocorrelation must be non negative. If it is zero then fog t dt has uniformly bounded 
L 2 norm. This implies that there exists a weakly converging subsequence to some L 2 function 
h such that E(/i) = 0. It is easy to check that such a function is smooth in the stable direction 
(just compare with the average on stable manifolds) and, for each smooth function ip, K(hL<p) = 
-E(ftp). Thus E(hL n tp) = -E(/L™-V) = (-1)"E(I™ -1 / (p), which implies L n h e L 2 , i.e. h is 
smooth along weak-stable leaves. Next, letting 0(q, v) = (q, —v), we have HL(ftp) = E(/i o Q • Lip), 
that is E((/i + ho Q)Lip) = for each smooth ip. In turns, this implies h = — h o a.s.. Indeed, 
given p S L 2 , if E(p) = and E(pLip) = for all smooth ip, then one can choose smooth p n that 
converges to p in L 2 , thus Lp n converges weakly to zero, but then there exist convex combinations 
p n of the {pm}m<n such that Lp n converges to zero strongly (since the weak closure of a convex 
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Next, consider the case r > 1. By Lemma 6.1 we have 

(B.2) T(r) = p(r, 1) = r" l p(l, r" 1 ) = -t^t" 1 , 1) = ^(t- 1 ). 

Thus r(r) = + 0(t~ 4 ). This readily implies the first part of the Lemma. 
Let us compute the derivative 

r'(r) 



dt m (LiV ■ L\V o g Tt ® g l ) 
oo /*rt 

dtt ds {E (ixV • L\V o 5 S ® 5 *) - E ■ L?F o 5 S ) } 
Jo 

+ / dt t I dsE {L{V ■ LfV o 5 s ) 
Jo Jo 

/ds dtt'E,(L 1 V-LlVog s ®g t )+ dt—E^iV-LlVogt) 

JQ Jt- 1 s ^ ' JO T 



oo c oo , /-oo u/r.i;. r ^ t/ 

ds 



o 



' dttE^LxV ■ LfV oid® g 1 ^) - j dt 



E (LiV • iiV o 5 * 



OO /-OO 



ds y df t y drE (LiV • if V" o g r ® 5 *) 

r r 00 / 9 ~ a E(iiF)E(F) -E(^i^) 

Jt- 1 s ^ 



OO 

1 



O I / dsr"^^"^ 



o 



= - r y dt t 2 E (l 2 ^ . r^y Q w 8 ^ + ( r 2) = . ^ + ^ r 2) 
On the other hand, differentiating (B.2) yields, for r large, 

r'(r) = -3T- 4 r(T- x ) - r-'^r'^- 1 ) = -3Ar- 4 + 0(t~ 5 ) 

which completes the proof of the Lemma. □ 

Remark B.2. Note that T(0) is not defined as the corresponding integral diverges. 
Nevertheless, we can set T(0) = A by continuity. 

Proof of Lemma 6.3. Note that, by Lemma 6.1, 

p(a, b) = b^piab- 1 , 1) = ft" 1 !^" 1 ). 

Hence the Lemma follows from Lemma B.l. □ 



set agrees with its strong closure) and, since L is a closed operator on L 2 , it follows that p is in 
the domain of L and Lp = 0. In addition, the ergodicity of the flow implies that the only L 2 , zero 
average, solution of Lp = is p = 0. Finally, since h is smooth along the weak-stable foliation 
and h o © is smooth along the unstable foliation, then h has a continuos version by the absolute 
continuity of the foliations and is smooth by [17], hence Lh = /. That is, if the autocorrelation 
is zero, then / is a smooth coboundary. At last, the claim follows since a smooth function of the 
coordinates only which is a coboundary must be identically zero, [10, Corollary 1.4]. Accordingly, 
J^^dt E (d qi V (<ji, q2)dq\ V(qi , g t (q2, ^2)) I must be strictly positive for positive measure 

set of qi otherwise, by the symmetry of the potential, the potential would be constant. 
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